tidystep - 分布

Code

library(gapminder)
library(tidyverse)
library(plotly)
library(scales)

Code

gapminder_data=gapminder 
glimpse(gapminder_data)

Rows: 1,704
Columns: 6
$ country   <fct> "Afghanistan", "Afghanistan", "Afghanistan", "Afghanistan", …
$ continent <fct> Asia, Asia, Asia, Asia, Asia, Asia, Asia, Asia, Asia, Asia, …
$ year      <int> 1952, 1957, 1962, 1967, 1972, 1977, 1982, 1987, 1992, 1997, …
$ lifeExp   <dbl> 28.801, 30.332, 31.997, 34.020, 36.088, 38.438, 39.854, 40.8…
$ pop       <int> 8425333, 9240934, 10267083, 11537966, 13079460, 14880372, 12…
$ gdpPercap <dbl> 779.4453, 820.8530, 853.1007, 836.1971, 739.9811, 786.1134, …

line

Code

data001=gapminder_data %>% group_by(year) %>% summarise(gdp=sum(gdpPercap)*sum(pop))
ggplot(data=data001 , aes(x=year, y=gdp)) +
  geom_point()+ geom_line()+scale_y_continuous(labels = comma)

point

Code

data001=gapminder_data %>% group_by(year,country) %>% summarise(gdp=sum(gdpPercap)*sum(pop))

ggplot(data=data001 , aes(x=year, y=gdp)) +
  geom_point()+scale_y_continuous(labels = comma)

boxplot

Code

data001=gapminder_data %>% group_by(year,country) %>% summarise(gdp=sum(gdpPercap)*sum(pop)) %>% filter(year>1990)

ggplot(data=data001 , aes(x=year, y=gdp,group=year)) +
  geom_boxplot()+geom_jitter()

quantile

Code

data001=gapminder_data %>% group_by(year,country) %>% summarise(gdp=sum(gdpPercap)*sum(pop)) %>% filter(year>1990)

data001  %>%  group_by(year) %>%
  summarize(quant10 = quantile(gdp, probs = 0.1)
,quant50 = quantile(gdp, probs = 0.5)
,quant90 = quantile(gdp, probs = 0.9)
)

# A tibble: 4 × 4
   year     quant10      quant50       quant90
  <int>       <dbl>        <dbl>         <dbl>
1  1992 3901409100. 34615235381. 435413157535.
2  1997 4166748953. 37473838315. 521107210532.
3  2002 4669205038. 41675762799. 605895961173.
4  2007 5918928216. 57869055458. 782717667549.

Code

gapminder_data=gapminder 
glimpse(gapminder_data)

Rows: 1,704
Columns: 6
$ country   <fct> "Afghanistan", "Afghanistan", "Afghanistan", "Afghanistan", …
$ continent <fct> Asia, Asia, Asia, Asia, Asia, Asia, Asia, Asia, Asia, Asia, …
$ year      <int> 1952, 1957, 1962, 1967, 1972, 1977, 1982, 1987, 1992, 1997, …
$ lifeExp   <dbl> 28.801, 30.332, 31.997, 34.020, 36.088, 38.438, 39.854, 40.8…
$ pop       <int> 8425333, 9240934, 10267083, 11537966, 13079460, 14880372, 12…
$ gdpPercap <dbl> 779.4453, 820.8530, 853.1007, 836.1971, 739.9811, 786.1134, …

Code

gapminder_data1997=gapminder_data %>% filter(year==1997)

Code

p <- ggplot(gapminder_data1997, aes(x=gdpPercap)) + 
  geom_histogram()
p

Code

p <- ggplot(gapminder_data, aes(x=gdpPercap,color=continent)) + 
  geom_histogram(fill="white", alpha=0.5, position="identity")
 
p

Code

p <- ggplot(gapminder_data, aes(x=gdpPercap,color=continent)) + 
  geom_histogram(fill="white", alpha=0.5, position="identity")+
  scale_x_continuous(trans='log10',labels = comma)
p

Reference